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' Abstract 

(N ' 

So how did the Galactic disk form and can the sequence of events ever be unravelled from the vast 
stellar inventory? This will require that some of the residual inhomogeneities from prehistory escaped the 
dissipative process at an early stage. Fossil hunting to date has concentrated mostly on the stellar halo, but 
<H ' a key source of information will be the thick disk. This is believed to be a 'snap frozen' relic which formed 

I ■ during or shortly after the last major epoch of dissipation, or it may have formed from infalling systems 

^•f) ' early in the life of the Galaxy. As part of the KAOS Galaxy Genesis project, we explore the early history of 

the halo and the thick disk by looking for discrete substructures, either due to infall or in situ star formation, 
t-H ■ through chemical tagging. This will require high signal to noise, echelle spectroscopy of up to a million stars 

k> throughout the disk. Our program has a short-term and a long-term goal. 

The short-term goal is to quantify the size and structure of the multi-dimensional chemical abundance 
space (C-space) for all major components of the Galaxy. We seek to establish how many axes in C-space 
are decoupled and have large intrinsic dispersions. A critical test of chemical tagging in the short term is 
that stellar streams in the halo, identified from detailed phase space information, are highly localized in 
C-space, or are confined to chemical tracks. These trajectories presuppose that stars form in a closed box 
\ through progressive enrichments of the gas, leading to stars dispersed along a narrow track in a complex 

chemical space. The long-term goal is to identify unique chemical signatures in the thick disk, originating 
' from different formation sites, for star clusters which have long since dispersed. This will require precise 

chemical abundances for heavy elements such that a star can be localized to a discrete point in C-space. If 
the star clusters originally formed outside the Galaxy in a bound infalling system, the stellar abundances 
' may fall along a chemical track, rather than a discrete point in C-space. 

C/2 . Keywords: Galaxy: kinematics and dynamics 
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1 Introduction turn. Many disk galaxies show a second fainter disk com- 

. £L ' . , , _ . , , fir ■ ponent with a longer scale heig ht (typically about 1 kpc) ; 

At the heart of a successful theory of gafaxy formation must ,,. . , ,, I„ , , ,. „ — rr~" 

, , ,. r i i- ■ ■ i this is known as the thi ck disk ( Schwarzkopf & Dcttmar , 

rN be a detailed physical understanding ol the dissipational \„ nnr l J77~i — — I — 7i — ~. — I n^T^T rTi TTTi TTt i 

, . , / . , , . 2000; Dalcanton & Bernstein . 2002). The Milky Way has 
processes which form spiral galaxies, the disk is the denn- .„ , , . . . ,. , „ I |-,„onl\ -j. r 

- , ? , a significant thick disk (Gilmorc & Kcid . 1983): its surface 

ing stellar component of disk galaxies, and understanding , . , , . , , „ r(w c "7T ,. , , c , . , , 

, r brightness is about 25% of the thm disk s surface bright- 

its formation is in our view the most important goal of , ., . .„ , , 

, r i r ness, and its stars are significantly more metal poor than 
galaxy formation theory. Although much of the mforma- 
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tion about the pre-disk state of the baryons has been lost 



the stars of the thin disk. 



The galactic thick disk is currently believed to arise 

m the dissipative process, some tracers are likely to remain. c , , , ,, ,. , , 

TTT , f i i ■ r i i from heating of the early stellar disk by one or more accre- 

What do we mean by the reconstruction of early galactic , . ... , ., , . . , , 

„ TTT , , ,. r i f lon events although other possible origins have been dis- 

history '. We seek a detailed physical understanding of the , 177^ — — 7T, „ T l ... 1 | lr , '^J 177^ — — 77; : — ^ — 

. * r .n ttt t, cussed (Gilmorc, Wvsc, & Kunkcn , 1989; Gilmorc. Wvse, & Jones , 

sequence of events which led to the Milky Way. Ideally, k „ nt i T , . , r , , , ... 

^ . r 1 ^ 1 11995)). It is composed of only old stars, with ages greater 

we would want to associate components of the Galaxy to , , . „ „ . , ,. . , , . r . , 

, than 10 Gyr, equivalent to lormmg at a redshitt 01 z > 1 

elements of the protocloud - the baryon reservoir which rrr 1 U„„A n ,1 ,1 , t • 1 1- 1 ■ 1 ■ n 

(Wvsc , 2000). hurthermore, the thick disk is chemically 

fueled the stars m the Galaxy. ,. , . , c ., ,,. , . , , „. .. JTTT - 

, distinct from the thm disk at low metalhcity (huhrmann . 

For many halo stars, and some outer bulge stars, this T7T7^\ U77T~- 7^ \ — r TT 1 I7TT^K ■ , 1 

, . , , . . . _ ' f998; FcltzmE, Bcnsbv, & Lundstrom , 2003), m the sense 

may be possible with phase space information. But for .. . r , . , , , , . . ... ljn , 

, , that a/he is enhanced relative to the thm disk, although 

much of the bulge and the disk, secular processes cause the , , , , . , . „. ., , ,. , 

11,1 1 r solar values are seen at higher metalhcity. the thick disk 

populations to become relaxed (i.e. the integrals of motion , rj , , . . n r , . 

v may be one of the most significant components for studying 
are partially randomized). In order to have any chance . r 1 r ,■ 1 •, > 

7 signatures of galaxy formation because it presents a snap 

of unravelling disk formation, we must explore chemical c , ,. ,,, . ,,, „ , , ,. , 

frozen refic of the state of the (heated) early disk, 
signatures m the stellar spectrum. T ,. . , . , r , ,. 

, ...... In the discussion which toflows, we make reference to 

The exponential thm disk, with a vertical scale height of , , . . , , r— n — — 1 — n — 1 ■ 1 

the KAOS project (see http://www.noao.edu/kaosl which 

about 300 pc, is the most conspicuous component m edge- , , ,,M , ., — ^-r? TT — 

, , proposes a highly multiplexed wide-fiefd multi-object spec- 

on disk galaxies. The thm disk is believed to be the end 



product of the quiescent dissipation of most of the baryons 
and contains almost all of the baryonic angular momen- 



trograph for the Gemini Observatory. Full specifications 
for the proposed instrument can be found at the above web 
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site. 



2 Unravelling a dissipative process 

In order to follow the sequence of events involved in dissipa- 
tion, we propose that the critical components which need 
to be re-assembled are the individual star clusters which 
formed at each stage. Since most stars are born in dense 
clusters, the formation and evolution of galaxies today must 
involve millions of discrete cluster events throughout their 
history. We would like to establish the evolving mass func- 
tion of star clusters, their chemical composition, formation 
and survival rate as a function of cosmic time. Galaxy-wide 
enrichment from the fall-out of nuclear winds or mergers 
would be evident in the fossil record of reconstructed star 
clusters, assuming these provide an unbiased sampling of 
cosmic time regardless of the star formation history. 

But how are we to reconstruct star clusters which have 
long since dispersed? It will be necessary to tag individual 
stars to their parent cloud through unique chemical signa- 
tures shared by these stars, assuming these exist. 

We now discuss our basic strategy fo r 'chemical tag- 
ging' l|Freeman fc Bland-Hawthorn 1 120021 . hereafter FBH). 
High resolution spectroscopy at high signal-to-noise ratio 
of many stellar types reveals an extraordinarily complex 
pattern of spectral lines. The spectral lines carry key in- 
formation on element abundances that make up the stellar 
atmosphere. Many of these elements cannot arise through 
normal stellar evolution, and therefore must reflect condi- 
tions in the progenitor cloud at the time of its formation. 



2.1 Primary requirements of chemical tagging. 

Our long-term goal is to chemically tag stars into coeval 
groups, i.e. to identify individual members of star clusters 
which have long since dispersed. For unique chemical sig- 
natures to exist, there are several key requirements (FBH): 

1. Most stars must be born in dense star clusters. 

2. Most dense star clusters must be chemically uniform 
in key elements. 

3. Key chemical elements must reflect the cloud com- 
position of the progenitor cloud. 

4. Key chemical elements must not be rigidly coupled 
(i.e. vary in lock step), and there must be sufficient abun- 
dance dispersion in key elements to allow for unique groups 
(reflecting unique sites of formation) to be readily identi- 
fied. 

5. There must be a contiguous spectral window which 
contains the necessary information on key elements for chem- 
ical tagging. 

We discuss each of these requirements, or conditions, 
in turn. Conditions 1, 2, 3 and 5 appear to be supported 
by observation. Condition 4 is the most uncertain largely 
because stellar abundance surveys to date target either too 
few stars or too few chemical elements. 

Condition 1. Most stars are born within rich clusters 



on Hyades bv lConti et al. I <ll965f) . But, until recently, very 
little was known about detailed heavy element abundance 
work on open clusters with rigorous membership estab- 
lished by reliable astrometry or accurate radial velocities 
(< 0.5 km/s), so as to minimiz e 'pollution' fro m stars not 
associated with the cluster (e.e. iQuillen L 12002'!. 

There have been several recent studies of open clus- 
ters which have established that open clusters are indeed 
highly u niform. In a detailed analysis of 55 F — K dwarfs in 
Hyades, iPaulson. Sneden. fc Cochran I (120031) find that the 
abundance variations in Si, Ti, Na, Mg, Ca and Zn with 
respect to Fe, and in [Fe/H], are within the measurement 
errors (<0.04 dex). 

There is remarkably little theoretical work which ad- 
dresses the chemical uniformity of star clusters. In theoret- 
ical work on star formation, chemical uniformity is almost 
always a prior assumption. Are we to expect chemical uni- 
formity among low mass stars in an open cluster? To date, 
thi s question has not bee n adequately answered. 

iMcKee fc Tan I (I2002T) propose that high- mass stars form 
in the cores of strongly self-gravitating and turbulent gas 
clouds. Two possible routes to chemical uniformity is that 
all stars form at the same instance from a chemically uni- 
form cloud, or that the low mass stars form outside of the 
core shortly after the supernovae have uniformly enriched 
the cloud. The precise sequence of events which give rise 
to open clusters is a topic of great int erest and heated de- 
bate in contemporary astrophysics (e.g. lStahler. Palla. Ho L 
l200d) . 

Condition 3. It is generally believed that r-process ele- 
ments (e.g. Sm, Eu, Gd, Tb, Dy, Ho) cannot be formed dur- 
ing quiescent stellar evolution. While some doubts remain, 
the most likely site for the r-process appear s to be Type II 
supern ova (SN II ) , as originally sugg e sted b v lBurbidee et alTI 
il957h (see also IWallerstein et al. L Il997ft . Therefore, r- 
process elements measured from stellar atmospheres reflect 
conditions in the progenitor cloud. The same is believed to 
be true for most of the a elements since these are produced 
in the hydrostatic burning phase of the pre-supernova star. 

In contrast, the s-process elements (e.g. Sr, Zr, Ba, Ce, 
La, Pb) are thought to arise from the He-burning phase 
of intermediate to low mass (AGB) stars (M<10 M©), al- 
though at the lowest metalli cities, trace amounts are likely 
to ari se from high mass stars feurris et al. 1l200f tlR.auscher et al. L 
120021) . 

Condition 4. During the past four decades, evidence 
has gradually accumulated for a large dispersion in metal 
abundances [X/Fe] (particularly n-capture el ements) in low 
metallicity stars relative to solar abundances (IWa llerstein et al. I 



| 1963 | : | Pagellll965USpite fc Spite Lll978l:lTruran Lll98ltlLuck fc Bond I 
198a:lciavton Ul988l:lGilrov et al. lll98alMcWirliam et al. I 
Il995l: iNorris. Rvan. fc Beers I 119961 : iBurris et al. L l2000l) . 

Elements like Sr, Ba and Eu show a 300-fold dispersion 



(see Fig. 21 of IWallerstein et al. L fl997l> . although [a/Fe] 
dispersions are typically an order of magnitude smaller. 

These observations have been used to constrain detailed 
supernova models, which in turn show how different yields 
of many hundreds to ma ny thousands of stars (IClarke. Bonnelte fefla$cfibfrffl«t .on of progenitor mass, progenitor metal- 



2000:; ICarpenter I l2000h . This essential fact is supported by 
many studies from optical, infrared, millimeter and radio 
surveys. 

Condition 2. The widely held view that open clus- 
ters are chemically uniform can be traced to classic work 



licity, mass cut (what gets ejected compared to what falls 
back towards the compact central object), and detonation 
details. These models help to explain the smaller disper- 
sions in a elements: the a yields are not dependent on the 
mass cut or details of the fallback/explosion mechanism, 
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Figure 1: An illustration of how successive chemi- 
cal enrichments by supernovae can lead to low mass 
stars with a wide scatter in chemical abundances 
which converge to a universal value. The model as- 
sumes a constant star formation rate in a closed box. 
The stars are formed according to a Salpeter mass 
function and the yields are t aken from models by 
iTsuiimoto fc Shigevama I l)l998|) . Only 1% of the stars 
produced in the simulation are shown after the initial 
burst; we show all the stars at t=0. 



which leads to a smaller dispersion at low metallicity. 

Is there evidence that some elements are decoupled? 
The element abundances [X/Fel show three mai n peaks at 
Z~26, Z~52, and Z~78 JSneden et aTL l2001allbl . There 
have been suggestions that th e r-process gives rise to ran - 
dom abundance patterns fe.g. iGorielv fc Arnould L 11996 h 
but this is not supported by new observations of a few 
metal poor stars. Heavy r-process elements around the 
second peak compared to the Sun appear to show a uni- 
versal sawtooth pattern, al beit with some internal vari- 
ance, in the range 56< Z<75 lISjMdeji^^a^JioOOyCavreU 
But lOtsuki. Mathews, fc Kaiino I 



2TO1I: iHill et, alTl 12 



( 2003) stress that this does not imply a universality of abun- 
dances out s ide th is range, and in fact this is not observed. 
IHill et al~l <l2002f) find that the third peak and actinide el- 
ements (Z>90) are decoupled from elements in the second 
peak. There may be a substantial number of suitable el- 
ements (10-20) which could define a sufficiently large pa- 
rameter space. 

The very large scatter of [n- capture/Fel in ultra-meta l 
poor stars seen in Figure 12 of IWallerstein et al. I Jl997t) 
means that n-capture element abundances in ultra-metal 
poor stars are products of one or very few prior nucleosyn- 
thesis events that occurred in the very early, poorly mixed 
galactic halo, a theme that has been developed by many 
authors (e.g. iGilrov et al. I Il988l: lAudouze fc Silkl Il995t 



McWilliam et al 
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erations of supernovae. The upper envelope is determined 
by enrichment from low mass supernovae, the lower enve- 
lope from high mass supernovae. From left to right, the 
black squares indicate the number of enrichment events, 
i.e. 1, 10, 100, 1000 events. Just ten enrichments from 
high mass supernovae are sufficient to enrich a cloud to 
[Fc/H]= —2. Note the rapid convergence in this simple 
model above [Fe/H]= —1. 

In the simulation, we consider only core collapse su- 
pernovae; subsequent enrichment by Type la supernovae 
causes the converging stream to dip down above [Fe/H] = 
— 1. It is not clear at what [Fe/H] the r-process elements 
become s wamped by the ubiquito us Fe-group and s-process 
elements. iTravaglio et al. 1 11199 91 suggest that the s-process 
does not become significant until |F e/H]^ —1 because of 
the need for pre-exi sting seed nuclei JSpite fc Spite Ul978t 
iTruran et akLl200ll) . although fpagel fc Tautvaisiene I (119*9% 
argue for some s-process production at [Fe/H]~ —2.5. In 
any event, Figure makes clear that, at sufficiently high 
resolution, both the Sun and Arcturus reveal a rich network 
of heavy elements, many of which arise from the r-process. 

An important point to appreci ate is that the substan - 
tially smaller scatter in [a/Fe] (see ICarretta et al. L 12002). 
compared to heavy element scatter, presumably argues that 
the scatter is unlikely to arise from chemical differentia- 
tion, i.e. where metal yields from a given supernova are 
distributed in different amounts, compared to Fe, to the 
surrounding ISM. In other words, the mixing is sufficiently 
uniform such that any given cell receives the same pro- 
portion of Eu/Fe from the supernova. This puts limits on 
the d egree to which asymmetric supernovae iMaeda et al. L 
2002) must differentially enrich the surrounding gas. 

Condition 5. Any practical experiment involving chem- 
ical signatures in millions of stars will require a wide-field, 
multi-fibre spectrograph operating at echelle resolutions. 
Since rayleigh scattering renders optical fibres almost use- 
less much below 4000 A, it is important to establish that 
there exist contiguous spectral windows which contain abun- 
dance information on a large number of chemical elements. 
Since it will be necessary to maximize the number of ob- 
ject fibres at the detector, this information will need to 
exist within a single narrow band (~ A/35) echelle order. 
The optimal parameters for a 'galaxy genesis' machine are 
determined below. 

Do there exist narrow band spectral windows, with chem- 
ical information on dozens of elements, which could be uti- 
lized by a fibre-fed echelle spectrograph? Until recently, 
this has not been possible to answer objectively. Figure |3 
shows such a window that was revealed by this study. An 
optimal window needs information on more than one ion- 
ization state for as many elements as possible in order to 
accurately determine the stellar surface gravity. Various 
ionic transitions in combination must also provide accurate 
information on the stars luminosity and surface tempera- 
ture. Some regions of the spectrum are dominated by tel- 
luric features; in cool stars, molecular bands complicate the 
spectrum. 

We co ncentrate our analysis on digitized atlases of the 



Tsuiimoto, Shigevama, & Yoshii 
The simulation in FigureQis illustrative of the expected 
abundance variation in successive generations of low mass 
stars in a closed box from gas enriched by successive gen- 



Tsuiimoto Ul998llArgagfti«tf pbl!Hl=0.0) and for the luminous giant Arcturu s ([Fe/H]= 



2000). 



Our 



-0.56) ( McWilli am. Rich, fc Smecker-Hane I _ 
solar line list derives from data supplied by P. Hall (NOAO) 
transcribed from the lMoore. Minnaert. fc Houtgastl <ll966l) 
solar spectrum (2935-8770 A). The atlas includes informa- 
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Figure 2: A demonstration of where most of the chemical information lies in a highly resolved spectrum of 
Arcturus (left) and the Sun (right). The black lines in the top plots give the total ion count in a spectral window 
AA = A/35. The magenta lines show what fraction of this count is made up of elements with two ionization 
states. The lower plots separate these into element types. Note the rapid rise of heavy elements to the blue in 
both stars; note also that the solar spectral information extends to bluer wavelengths. 



Figure 3: 17 A spectral window (4047-4064 A) of 
the Sun obtained with the Fourier Transform Spec- 
trograph at the National Solar Observatory reveal- 
ing detailed chemistry on Fe, Cr, Ti, V, Co, Mg, 
Mn, Nd, Cu, Ce, Sc, Gd, Zr, Dy. The lower values 
are wavenumbers (cm -1 ). The data are taken from 
http: / /www. nso.edu/ 



tion on 25,000 absorption lines covering 63 elements, for 
which 28 are detected in two ionization states. The Arc- 
tu rus line list comes from the digi t ized a tlas (3570-7405A) 
of IWallace. Hinkle. fe Livingston! Jl998h . The line list in- 
cludes 7700 atomic transition for 42 elements, for which 18 
are detected in two ionization states. Both lists are domi- 
nated by Fe I (45%) and light a elements (30%). 

Figure |5] shows the number of unique ions detectable 
as a function of wavelength for the Sun and for Arcturus. 
At each wavelength, the accessible window is AA where 
A/AA =35, such that at A — 3500 A, the unique ion count 
is made over a 100 A. The linear dependence of the window 
on wavelength is expected for spectral coverage in a single 
echelle order limited by the detector size. Note the rising 
count towards with decreasing wavelength down to 4000 A, 
both in terms of the number of unique ions, and elements 
with two ionization states. 

In the lower panels of Figure [5] we demonstrate that 
this rise is largely due to a rising fraction of neutron capture 
elements; the rise in iron peak elements is not as dramatic. 
It is noteworthy that we obtain similar results for a solar 
metallicity main sequence star and a metal weak giant. Of 
particular importance to the KAOS project, note that there 
is an optimal window close to 4100 A accessible to optical 
fibres. 

In Figure 2] we show the effect of resolving power on 
the fraction of lines which can be resolved as a function 
of wavelength. The optimal window at 4100 A requires a 
spectral resolving power of R = 40,000 in order to resolve 
80 % of the lines. As the resolving power decreases, the loss 
of information on heavy elements is particularly dramatic. 

We conclude that the KAOS Galaxy Genesis project 
should target 4100 A in a 100 A window at a resolving 
power of R = 40,000. The required sampling 0.05 A per 
pixel or 2000 pixels in the dispersion direction, well matched 
to state of the art CCDs. 
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Figure 4: The effect of resolving power on the fraction of spectral lines which can be resolved as a function of 
wavelength. In order to resolve 80 % of lines at 4000 A, we need R = 40,000; this falls to R = 10,000 at 8000 
A. Once again, both stars reveal the same general trends; note also that the solar spectral information extends 
to bluer wavelengths. 



2.2 Candidates for chemical tagging 

Chemical tagging is not possible for all stars. In hot stars 
and young star clusters, our ability to measure abundances 
is reduced by the stellar rotation and lack of transitions for 
many ions in the optical. In very cool stars, very little light 
emerges in the optical or near-IR due to complex molec- 
ular and dust opacities. The ideal candidates are evolved 
FGK stars that are intrinsically bright, FGK subgiants and 
dwarfs. These populations probably account for 10% of all 
stars. 

Giants can be observed at R = 40,000 over the full 
Gaiasphere, i.e. the 20 kpc diameter sphere centred on 
the Sun which will be surveyed in great detail by the GAIA 
satellite. Dwarf stars will only be observable within 1-2 kpc. 
While this is only a small fraction of the available volume, 
the local volume may include a representative sample of all 
old disk stars, regardless of their point of origin. 

For many stars, certain spectral regions are dominated 
by thermally broadened profiles, particularly for light ele- 
ments. In general, heavier atoms show narrow profiles ap- 
propriate to their mass, although odd numbered atomic iso- 
topes are susceptible to hyperfine splitting (due to the non- 
zero nuclear magnetic moment) which produces a broad- 
ened asymmetric line profile. 

Giants have deep, low density atmospheres that pro- 
duce strong low-ionization absorption lines compared to 
higher gravity atmospheres. Even in the presence of sig- 
nificant line blending, with sufficient signal, it should be 
possible to derive abundance information by comparing the 
fine structure information with accurate stellar synthesis 
models. Detailed abundances of large numbers of F and G 
subgiants would be particularly useful, if it becomes possi- 
ble to make such studies, because direct relative ages can 
be derived for these stars from their observed luminosities. 



3 Short-term goal size and structure in a 
multi-dimensional C-space 

An intriguing prospect is that reconstructed star clusters 
can be placed into an evolutionary sequence, i.e. a family 
tree, based on their chemical signatures. Let us suppose 
that a star cluster has accurate chemical abundances de- 
termined for a large number n of elements (including iso- 
topes). This gives it a unique location in an n-dimensional 
space compared to m other star clusters within that space. 
We write the chemical abundance space as C(Fe/H, Xi/Fe, 
X2/Fe, ...) where Xi, X2 ... are the independent chemical 
elements that define the space (i.e. elements whose abun- 
dances are not rigidly coupled to other elements). 

Our simple picture assumes that a cloud forms with a 
unique chemical signature, or that shortly after the cloud 
collapses, one or two massive SN II enrich the cloud with 
unique yields which add to the existing chemical signature. 
The low-mass population forms with this unique chemical 
signature. If the star-formation efficiency is high (>30%), 
the star group stays bound although the remaining gas is 
blown away. If the star-formation efficiency is low, the star 
cluster disperses along with the gas. In a closed box model, 
the dispersed gas reforms a cloud at a later stage. 

In the closed box model, each successive generation of 
supernovae produce stellar populations with progressive en- 
richments . These will lie along a t r ajecto ry in C-space (see 
Fig. 8 of iKarlsson fc Gustafsson 1 1200 ll) . The overall dis- 
tribution of the trajectories will be affected by fundamental 
processes like the star formation efficiency, the star forma- 
tion timescale, the mixing efficiency, the mixing timescale, 
and the satellite galaxy infall rate. 

A critical test of chemical tagging is that stellar streams 
in the halo, identified from detailed phase space informa- 
tion, are highly localized in C-space, or are confined to 
chemical tracks this is a key short term goal. There may 
already be evidence for accreted halo stars from their dis- 
tinct chemical signatures (e.g. suppressed Mg) since those 
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with estimated orbital p aramet e rs are found to have large 
apogalacticon distances dKing L ll997L iNissen fc Schuster 1 
Il997t ICarnev et aTl Il997l : ICarretta et al. I l2002t) . This 
may suggest that these stars originated in lower mass stel- 
lar systems with very different chemical histories from that 
of the Galaxy. 

Our team (De Silva, Bland-Hawthorn, Freeman, Bessell, 
Asplund) is currently engaged in large surveys to establish 
the chemical homogeneity of open clusters as this forms 
a central tenet of chemical tagging. We propose to ex- 
tend this survey to large star-forming complexes which in- 
clude several stellar associations within one or two super- 
associations. This is the sample volume we need to tackle 
to establish the existence of chemical tracks. Presently, it is 
difficult to identify a proper sample until phase space infor- 
mation becomes available from GAIA, although Hipparcos 
has pro vided detailed candidacy for a few star formin g com- 
plexes dde Bruiine. Hooeerwerf. fc de Zeeuw L 1200 if) . 

As we approach solar levels of metallicity in [Fe/H], the 
vast number of trajector i es wil l tend to converge (see Fig. 

21 of IWallerstein et al. L fl997h . By [Fe/H] 2.5, AGB 

stars will have substantially raised the s-process element 
abundances; by [Fe/H]~ —1, Type la supernovae will have 
raised the Fe-group abundances. Star clusters that appear 
to originate at the same location in this C-space may simply 
reflect a common formation site, i.e. the resolution limit we 
can expect to achieve in configuration space. 

Even with a well established family tree based on chem- 
ical trajectories in the chemical C-space, this information 
may not give a clear indication of the original location 
within the protocloud or Galactic component. This will 
come in the future from realistic baryon dissipation mod- 
els. Forward evolution of any proposed model must be able 
to produce the observed chemical tree. 

However, the C-space will provide a vast amount of in- 
formation on chemical evolution history. It should be pos- 
sible to detect th e evolution of the cluster mass function 
with cosmic time (lKroupal i2002). the epoch of a starburst 
phase and /or as sociated mass ejection of metals to the hal o 
dRenzini L I2002T) . and/or satellite infall dNoguchi L ll99Sl> . 
The chemical tracks could conceivably be punctuated by 
discontinuities du e to dramatic events like galactic infall or 
large-scale winds dChiappini et al.1ll999Tl . 

As we go back in time to the formation of the disk, 
we approach the chemical state laid down by population 
III stars. The rarity of stars below [Fe/H]~ —5 suggests 
that the protocl oud was initial l y enri ched by the first gen- 
eration of stars dArgast et al. L 120001) or maybe that stars 
moving through the ISM today have a minimum thresh- 
old metallicity due to Bondi-Hoyle accretion. However, 
the apparent absence of any remnants of population III 
remains a puzzle: its stars may have had a top-heavy ini- 
tial mass function, or have dispersed into the intra-group 
medium of the Local Group. If one could unravel the 
abundances of heavy elements at the time of disk forma- 
tion, this would greatly improve the precision of nucleo- 
cosmochronology. Important information is beginning to 
emerge from ech elle observations of damped Lya systems 
at high redshift dPettini et al. 112003ft . 



3.1 The size of C-space 

The size of n is unlikely to exceed about 60 for the fore- 
seeable future. ISneden et al. I <|2003T) present exquisite data 
for the metal-poor star CS 22892-052, where abundance es- 
timates are obtained for a total of 54 elements, w ith tight 
upper limits on 3 other elements. iHill et al~l i2002fl present 
exquisite data for the metal-poor star CS 31082-001, where 
abundance estimates are obtained for a total of 44 elements , 
almost half the entire periodic table (see also lCavrel U200ll 
ISneden et al. llioOlal lbl). The a and r- process elements, and 
maybe a few canonical s-process elements at low [Fe/H], 
provide information on the cloud abundances prior to star 
formatio n, although combinations of these are likely to be 
coupled dHeger fc Wooslev U200 lllSneden et al. Ll2001allb|) . 
There are more than two dozen r-process elements that have 
been clearly identified in st ellar spectra dWallerstein et al~l 
ll997l:ISneden et al. LEoo3> . 

The size of m is likely to be exceedingly large for the 
thin disk where most of the baryons reside. For a rough es- 
timate, we take the age of the disk to be 10 Gyr. If there is 
a unique SN II enrichment event every 100 years, we expect 
of order 10 s formation sites. Typically, a SN II event sweeps 
up a constant mass of 5xl0 4 dRyan. Norris. fc Beers Ul99rj 
IShigevama fc TsuiimotoT I1998T) . Simple chemical evolu- 
tion models indicate that this must be of the ri ght order to 
expla in the metallicity dispersion at low [Fe/H] dArgast et al. I 
2000). Roughly speaking, there have been 10 3 generations 
of clouds since the disk formed, with about 10 5 clouds in 
each star-forming generation, such tha t cloud formation 
and d ispersal cycles on a 10 7 yr timescale dElmegreen et al~l 
120001) . 

Whereas the total number of star clusters over the life- 
time of the thin disk is very large , the size of m for the 
stella r halo l|Harding et al. L 1200 J) , and maybe the thick 
disk dKroupaL " l2002r i. is likely to be significantly smaller. 
Our primary interest is the oldest star clusters. Recon- 
structing star clusters within the thick disk is a particularly 
interesting prospect since the disk is likely to have formed 
within 1-1.5 Gyr of the main epoch of baryon dissipation 
dProchaska et al. Ll200oT) . 

3.2 How many unique chemical signatures make 
up the halo or the thick disk? 

There are different approaches to answering this question: 

(i) carry out a cold dark matter simulation (with hydro) 
and trace stellar evolution and feedback within each of the 
building blocks which make up a galaxy; (ii) carry out a 
Press-Schechter simulation based on statistical aggregation 
with prescriptive physical laws; (iii) undertake an order of 
magnitude calculation. 

Our basic view is that modern CDM simulations in (i) 
are not able to carry enough chemical information in order 
to do this properly. In its place, the numerical approach in 

(ii) is currently under investigation by us: our initial con- 
clusion is that the expected number of signatures is very 
poorly defined due to uncertainties in supernova yields, and 
the sequence of events involved in star formation and recy- 
cling of gas. In place of detailed simulations, we provide an 
order of magnitude estimate. 

There are at least three possible scenarios for the thick 
disk, all of which can be tested by chemical tagging. (1) 
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The thick disk is thought to be a snap-frozen relic of the 
early disk, heated vertically by the infall of an intermediate 
mass satellite. This is a particularly interesting prospect 
since chemical tagging would provide clues on t he forma- 
tion o f the first star clusters in the early disk. (2l lKroupal 
(2002) has suggested that the thick disk arose from a pop- 
ulation of su per-star clusters ( ~ 10 6 Mq) which became 
unbound. (3) iBekki fc Chiba I j200Cf) have shown how the 
thick disk could arise fro m the infall o f a num ber of Local 
Group satellites; see also lAbadi et al. I J2003ft . Their mod- 
els could explain why thick disk stars are more metal rich 
than stars in dwarf galaxies or other likely building blocks, 
because the minor merger would have triggered star forma- 
tion and self-enrichment. 

For a thick disk made up of super-star clusters, we 
would need to detect roughly 10 4 unique chemical signa- 
tures, requiring a survey of 10 stars in order to detect a 
unique chemical signature at the 3a level. For a thick disk 
made up of the heated early disk, we would need to detect 
of order 10 7 unique signatures, requiring a survey of 10 s 
stars. For a thick disk comprising roughly ten Local Group 
satellites, the number of unique chemical signatures will 
depend on the details of the star formation history within 
these satellites prior to infall. But very roughly, we would 
need to detect 10 6 unique signatures, requiring a survey of 
10 7 stars. 

In the next section, in a more detailed statistical analy- 
sis, we show that a survey of 10 6 stars may be sufficient to 
identify the chemical signatures of the early thick disk and 
the stellar halo. 

3.3 Searching for progenitor formation sites — 
how many stars do we need? 

Our simulations show that we need to identify roughly 10% 
of the original formation sites in order to properly sam- 
ple cosmic time for a wide range of possible star formation 
histories (stochastic, exponential). In the table below, we 
show the size of survey needed to detect a certain fraction 
of formation sites where the detection requires 1, 2, 3, 5, 10 
or 30 stars from each site (T able |T) l. 

Thick disk. If we adopt iKroupa I i2002f) model for the 
thick disk, we would need to survey 60,000 thick disk stars 
in order to establish 3<r chemical signatures for 10% of the 
original formation sites. As we show below, this requires 
500—1000 KAOS pointings in order to survey such a large 
number of thick disk stars in a limited magnitude range. 

Stellar halo. The baryon mass of the stellar halo is in 
the range 10 8 to 10 9 Mq depending on the details of how 
we distinguish between halo and bulge stars, and whether 
or not we include recently identified halo objects (e.g. Sgr 
dwarf). For illustration, we take the lower bound here since 
we can scale trivially to the higher bound. 

Early universe simulations suggest that the first star 
clusters (e.g. protoglobulars, first dwarfs) mus t contain 
rough ly a Jeans mass, or about 10 5 Mq in baryons jAbel 
2002). Since a single supernova can enrich a cloud of 10 5 
Mq, let us suppose that the Galactic halo arises from 1000 
progenitor sites. Therefore, we would need to survey 500 
halo stars to be sure of seeing the chemical signatures of 
more than one star from 10 % of the original formation 
sites. With sufficiently accurate abundance information, a 



survey of 6000 halo stars would give a discrete chemical 
signatures (3<r) for 10 % of the original sites. 

If the original halo formation sites saw several genera- 
tions of star formation before infall, we should expect to 
see chemical tracks in C-space. If we need 30 stars from a 
single formation site to establish a chemical track, we need 
to survey roughly 25,000 halo stars in order to establish 
tracks for 10% of the progenitor sites. This amounts to 
500-1000 KAOS pointings which can be observed simulta- 
neously with the thick disk pointings (see below). 

If chemical tracks exist, it is fundamentally important 
to identify them because the starting point identifies a unique 
set of chemical signatures from a first-generation super- 
nova. As mentioned above, there are a host of fundamen- 
tal uncertainties about the enrichment details from super- 
novae, which severely limits our ability to perform real- 
istic CDM+hydro simulations. But there is likely to be 
substantial progress in this arena Ce.g. lGibson et a! U20o3 
iBrook et all I2004T) . 

Note that the chemical tagging of the halo will be as- 
sisted by phase space information supplied by the KAOS 
survey and the GAIA space mission. It is presently unclear 
whether phase space information for the thick disk will pro- 
vide supporting evidence for discrete groups identified by 
their unique chemical signatures. 

In summary, a survey of 25,000 halo stars and 60,000 
thick disk stars may be sufficient to detect a useful number 
of progenitor formation sites; these observations will share 
500-1000 KAOS pointings (see below). 

3.4 How many chemical signatures are we ever 
likely to detect? 

The data to answer this question do not exist: stellar sur- 
veys to date target either too few stars or too few chemical 
elements. The answer depends in part on the strategy: we 
could attempt accurate measurements (<0.05 dex) on a few 
decoupled lines, or try to distinguish, say, a high and a low 
state in many elements. An error of 0.02 dex is equivalent 
to 1 mA measurement accuracy in equivalent width. 

Most physical processes are intrinsically noisy, so it 
would be surprising if there did not exist at least a small 
amount of decoupling in elements which appear to corre- 
late within the measurement errors over a wide dynamic 
range in metallicity. The rapid statistica l convergence with 
increa sing [Fe/H] seen in Figure 21 of IWallerstein et al. I 
lll997l) is assumed to arise from the increasing homogene- 
ity of the ISM with cosmic time. But some intri nsic scatter 
could be sustained by gaseous infall (however see lPaeel fc Tautvaisiene"! 
I1997D . 

In Figure [5] there are more than 30 unique elements 
observable in our optimal window. If we could distinguish 
a high and a low state for each of these, we could in princi- 
ple distinguish a billion unique chemical signatures. Either 
strategy would require an accurate differential abundance 
•aansflyMBr(fiflliv;irdsson et al. 1 119931 iProchaska et al. 1 l2000t 
iPaulson. Sneden. fc Cochran U2003I) . 

But we stress that it may not be necessary to measure as 
many as 30 elements if some can be found which are highly 
decoupl ed and exhibit large relative dispersions from star 
to star. iBurris et al. I (l2000f) demonstrate o ne such element 
pair, i.e. [Ba/Fe] vs [Sr/Fe], confirmed bv ICarretta et al~l 
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Table 1: The number of stars needed to detect progenitor sites (normalized to the number of sites). The rows 
indicate the fraction of sites which are detected in 1, 2, 3, 5, 10, 30 stars. For example, if there are 1000 possible 
sites, and we want to detect 5 stars in 10% of these, we need to randomly sample about 2400 stars from the 
complete sample. 





10% 


20% 


30% 


40% 


50% 


60% 


70% 


80% 


90% 


1 


0.11 


0.22 


0.36 


0.50 


0.68 


0.90 


1.19 


1.60 


2.32 


2 


0.53 


0.81 


1.08 


1.36 


1.67 


1.99 


2.41 


2.97 


3.87 


3 


1.10 


1.53 


1.92 


2.27 


2.67 


3.10 


3.62 


4.26 


5.29 


5 


2.43 


3.09 


3.61 


4.13 


4.65 


5.21 


5.88 


6.73 


8.08 


10 


6.22 


7.30 


8.16 


8.90 


9.65 


10.46 


11.41 


12.60 


14.21 


30 


23.20 


25.30 


26.92 


28.33 


29.64 


31.05 


32.69 


34.56 


37.21 



(2002) who also suggest [Mn/Fe] vs [Cr/Fe] may have ap- 
preciable dispersion. It is part of our short-term goal to es- 
tablish which axes in C-space have the largest dispersions, 
and which axes are likely to be most decoupled. 

4 Long-term goal reconstructing ancient star 
groups from unique chemical signatures 

The abundance dispersion in a and heavy elements pro- 
vides a route forward for tagging groups of stars to common 
sites of formation. With sufficiently detailed spectral line 
information, it is feasible that the 'chemical tagging' will 
allow temporal sequencing of a large fraction of stars in a 
manner analogous to building a family tree through DNA 
sequencing. 

Consider the (extraordinary) possibility that we could 
put many coeval star groups back together over the entire 
age of the Galaxy. This would provide an accurate age 
for the star groups either through the color-magnitude di- 
agram, or through association with those stars within each 
group that have [n-capture/Fe] 2> 0, and can therefore be 
radioactively dated. This would provide key information 
on the chemical evolution history for each of the main com- 
ponents of the Galaxy. 

There is no known age-metallicity relation that operates 
over a useful dynamic range in age and/or metallicity. (This 
effect is only seen in a small subset of hot metal- rich stars). 
Such a relation would require the metals to be well mixed 
over large volumes of the ISM. For the foreseeable future, 
it seems that only a small fraction of stars can be dated 
directly (FBH). 

Ideally, we would like to tag a large sample of represen- 
tative stars with a precise time and a precise site of forma- 
tion. Can we identify the formation site? The kinematic 
signatures will identify which component of the Galaxy the 
reconstructed star group belongs, but not specifically where 
in the Galactic component (e.g. radius) the star group came 
into existence. For stars in the thin disk and bulge, the stel- 
lar kinematics will have been much affected by the bar and 
spiral waves; it will no longer be possible to estimate their 
birthplace from their kinematics. Our expectation is that 
the derived family tree will severely restrict the possible sce- 
narios involved in the dissipation process. In this respect, a 
sufficiently detailed model may be able to locate each star 
group within the simulated time sequence. 



Our ability to detect structure in C-space depends on 
how precisely we can measure abundance differences be- 
tween stars. It may be possible to construct a large database 
of differential abundances from echelle spectra, with a pre- 
cision of 0.05 dex or better; differential abundances are pre- 
ferred here to reduce the effects of systematic error. 

4.1 The observing programme 

We seek to target old stars in the thick disk and halo. In 
order to target a large enough volume of the thick disk and 
halo, we need reach to an apparent mag V = 17 (equiv- 
alent to I ~ 16.5); see Table H We propose to observe 
stars in a narrow magnitude range, 15 < V < 17, in order 
to reduce the effects of scatte red light in the instrument 

JWvse fc GilmoreUl99lll995h 

If we adopt Gilmore's model llGilmore. Reid. fc Hewitt 1 

Il985l) of the Galaxy (see also iRobin et al. I l2003fi . within 
our magnitude range towards the Galactic poles, there are 
about 230 stars deg -2 with the following breakdown: 140 
thin disk main sequence (MS), 45 thick disk MS, 10 halo 
MS, 15 evolved (subgiant or red giant) thick disk, 15 evolved 
halo, 5 halo horizontal branch. Along a cardinal sight line of 
(1=90°, b=30°), these numbers increase to 820 stars deg -2 , 
with the following breakdown: 660 thin disk main sequence 
(MS), 80 thick disk MS, 10 halo MS, 4 evolved thin disk, 43 
evolved (subgiant or red giant) thick disk, 15 evolved halo, 
5 halo horizontal branch. This amounts to 60 thick disk 
stars towards the poles, and 120 stars at the lower latitude; 
we observe about 30 halo stars along either sight line. 

A reasoned SNR calculation for a high-resolution KAOS 
design (see web calculator at http://www.noao.edu/kaos) 
shows that we should be able to achieve SNR~100 per pixel 
at V = 17 in a 4 hr exposure (R~20,000) rising to SNR-300 
at V = 15. Here we assume a 0.5" seeing observation of a 
3000K blackbody source at 5000 A in full moon. A me- 
dian SNR~150 for the survey will be sufficient to identify 
unique chemical tracks and chemical signatures. This level 
of performance is encouraging for R~40,000 performance if 
this could be factored into the KAOS design. 

For sanity, we have checked the predicted performance 
against a UVES commissioning observation (single slit, not 
fibre fed) which claims SNR~100 at 4000 A for I~15 star 
at R~25,000 in 1 hour exposure. This is a factor of two 
better than KAOS as expected since KAOS is fibre fed. 

The actual spectroscopic resolution we need to achieve 
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Table 2: Distance limit in log(parsecs) for different stars as a function of apparent V magnitude: metal poor 
giants (MPG), metal rich giants (MRG), clump giants (CG), blue horizontal branch halo (BHB), and main 
sequence dwarfs. The second column is the absolute V magnitude of the star. Brackets help to delineate the 
transition between 1 - 10 - 100 kpc. Note that the Solar Circle provides an extra 8 kpc in radial extent such 
that surveys which reach the Galactic Center also reach the outer disk. 





V 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


MPG 


-2.0 


(4.0) 


4.2 


4.4 


4.6 


4.8 


(5.0) 


5.2 


5.4 




-1.5 


3.9 


4.1 


4.3 


4.5 


4.7 


4.9 


5.1 


5.3 


MRG 


-1.0 


3.8 


(4.0) 


4.2 


4.4 


4.6 


4.8 


(5.0) 


5.2 




-0.5 


3.7 


3.9 


4.1 


4.3 


4.5 


4.7 


4.9 


5.1 




0.0 


3.6 


3.8 


(4.0) 


4.2 


4.4 


4.6 


4.8 


(5.0) 


CG/BHB 


0.5 


3.5 


3.7 


3.9 


4.1 


4.3 


4.5 


4.7 


4.9 




1.0 


3.4 


3.6 


3.8 


(4.0) 


4.2 


4.4 


4.6 


4.8 




1.5 


3.3 


3.5 


3.7 


3.9 


4.1 


4.3 


4.5 


4.7 


A 


2.0 


3.2 


3.4 


3.6 


3.8 


(4.0) 


4.2 


4.4 


4.6 




2.5 


3.1 


3.3 


3.5 


3.7 


3.9 


4.1 


4.3 


4.5 




3.0 


(3.0) 


3.2 


3.4 


3.6 


3.8 


(4.0) 


4.2 


4.4 


F 


3.5 


2.9 


3.1 


3.3 


3.5 


3.7 


3.9 


4.1 


4.3 




4.0 


2.8 


(3.0) 


3.2 


3.4 


3.6 


3.8 


(4.0) 


4.2 




4.5 


2.7 


2.9 


3.1 


3.3 


3.5 


3.7 


3.9 


4.1 


G 


5.0 


2.6 


2.8 


(3.0) 


3.2 


3.4 


3.6 


3.8 


(4.0) 




5.5 


2.5 


2.7 


2.9 


3.1 


3.3 


3.5 


3.7 


3.9 




6.0 


2.4 


2.6 


2.8 


(3.0) 


3.2 


3.4 


3.6 


3.8 




G.5 


2.3 


2.5 


2.7 


2.9 


3.1 


3.3 


3.5 


3.7 


K 


7.0 


2.2 


2.4 


2.6 


2.8 


(3.0) 


3.2 


3.4 


3.6 




7.5 


2.1 


2.3 


2.5 


2.7 


2.9 


3.1 


3.3 


3.5 



depends in part on the spectrophotometric stability, both 
in terms of an accurate zero level to each spectrum, and 
an accurately calibrated instrument response. If both zero 
level and instrument response can be calibrated accurately 
and consistently, it may be possible to conduct the KAOS 
survey at R = 20,000. However, if, say, the zero level can 
be maintained, but not the instrument response correction, 
it is safer to work at R = 40,000, and fit individual lines 
in order to measure equivalent widths. But the advantage 
of removing the instrument response accurately is that it 
assists in determining the true continuum level by matching 
to atmospheric models. If more than one order is observed, 
instrument response correction is also important for joining 
orders together. 

In summary, for the benchmark KAOS one degree field, 
we detect 60,000 thick disk stars in 500-1000 KAOS fields, 
and 25,000 halo stars in about 1000 KAOS fields. Ideally, 
we would opt for R~40,000 but it may be possible to work 
at R~20,000 with sufficient instrument stability. 



5 Related surveys and the way ahead for near- 
field cosmology 

Detailed high resolution abundance studies of large samples 
of galactic stars will be crucial to the success of the Galaxy 
Genesis project. The KAOS survey is one of three massive 
stellar surveys. 

The first is the RAVE survey (see http: / /www. iap.de/R7 
which will obtain radial velocities and abundances for 50 



million stars, i. e. all stars down to V=16, over the pe- 
riod 2005-2010 JSteinmetz I l2002h . This requires a new 
generation 'Echidna' spectrograph mounted on a Schmidt 
telescope (1.2m aperture), with 2250 fibres working at a 
spectroscopic resolution of R ~ 8000. 

By far the most ambitious proposal is the GAIA satel- 
lite mission (see http://www.rssd.esa.int/GAIA) over the 
period 2012-2018 which is expect ed to launch dur i ng th e 
time frame of the KAOS survey iPerrvman et al. 1 12001 ). 
This will provide photometric, spectroscopic and phase space 
information for a billion stars in the Local Group. The 
GAIA mission is expected to obtain R ~ 10,000 resolution 
spectroscopy for up to 100 million stars or more. 

It is important to note several points about the GAIA 
and RAVE surveys: these are based on snapshot surveys 
with a lm aperture telescope, targeting the 8700A region 
at low resolving power. Their choice of R = 8000 (RAVE) 
and R = 10,000 (GAIA) are supported by Figure El As 
we have seen, this region carries limited abundance infor- 
mation. Both missions are primarily optimized for radial 
velocity studies. The KAOS echelle survey will be the first 
million star survey optimized for chemical abundance work. 

The KAOS survey requires a highly multiplexed, high 
resolution spectrograph which will be expensive and tech- 
nically challenging, but we believe this must be tackled if 
we are to ever unravel the formation of the Galaxy. Could 
some of the residual inhomogeneities from prehistory have 
escaped the dissipative process at an early stage? We may 
r &k, know the answer to this question with absolute cer- 
tainty for many years. But it is an intriguing thought that 
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one day we may be able to identify the Solar Family, i.e. 
the hundreds or even thousands of stars throughout the 
Gaiasphere that were born within the same cloud as the 
Sun. 
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